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Abstract 



Elastic and proton-dissociative photoproduction of p , <f) and J/ip vector mesons (jp —>■ Vp, 

jp — > VN, respectively) have been measured in e + p interactions at HERA up to — t = 3 GeV 2 , 

where t is the four-momentum transfer squared at the photon-vector meson vertex. The analysis 

i-J , is based on a data sample in which photoproduction reactions were tagged by detection of 

the scattered positron in a special-purpose calorimeter. This limits the photon virtuality, Q 2 , 

to values less than 0.01 GeV 2 , and selects a jp average center-of-mass energy of (W) = 94 

OS ' GeV. Results for the differential cross sections, da/dt, for p°, <\> and J/ip mesons are presented 

and compared to the results of recent QCD calculations. Results are also presented for the 

^ ' ^-dependence of the pion-pair invariant-mass distribution in the p° mass region and of the spin- 

1} , density matrix elements determined from the decay-angle distributions. The Pomeron trajectory 

has been derived from measurements of the W dependence of the elastic differential cross sections 

D , da/dt for both p° and <f> mesons. 
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1 Introduction 

The study of exclusive diffractive ep reactions at HERA has shown that when Q 2 (photon virtuality) 
or the mass scale involved is large, the cross section increases with energy faster than expected for 
soft processes Q ^] . The rise is consistent with predictions from models based on perturbative 
QCD (pQCD) in which Q 2 and mass are used as the perturbative scale. It is expected that, for 
diffractive vector-meson photoproduction, the four-momentum transfer squared, t, between the 
photon and the final-state vector meson may also serve as a hard scale, provided —t is large H. In 
the present paper, this hypothesis is studied by measuring diffractive vector-meson photoproduction 
as a function of t and by confronting the data with the predictions of models based on pQCD, which 
should be applicable in the presence of a hard scale. In addition, the data at lower —t are compared 
to predictions of models expected to be valid for soft processes. This provides a means to study 
the transition between the soft, non-perturbative, and the hard, perturbative, regimes of QCD 0]. 
A detailed study of the onset of the pQCD regime should give important insight into the structure 
of strong interactions at hard scales as well as improve our understanding of soft phenomena in 
QCD. 

Vector mesons can be diffractively photoproduced via two processes. In one of them, the target 
proton remains intact and the reaction is called exclusive (or elastic), 

7P — > Vp. (1) 

In the other process, the proton dissociates into a higher mass nucleonic state N and the reaction 
is called proton-dissociative, 

lp _> VN. (2) 

Reaction (|l|) is called elastic in the framework of the vector-dominance model (VDM) ||, in which 
the photon fluctuates into a virtual vector meson which in turn scatters elastically from the target 
proton ||. This reaction has been studied over a wide range of 7p center-of-mass energies W < 
200 GeV |l], q. The proton-dissociative reaction (p|) has been studied at low energies and 
information at high W has been obtained recently J8|. 

In this study, reactions (||) and (g) are investigated with the ZEUS detector at HERA by measuring 
the processes ep — ► eVp and ep — ► eVN, where V = p°,4> or J/^i an d N is a system with 
mass < 7 GeV. The scattered positron was detected in an electromagnetic calorimeter close to the 
beamline at a distance of 44 m from the interaction point in the direction of the outgoing positron. 
This ensured that the virtuality of the exchanged photon is very small (Q 2 < 0.01 GeV 2 ) and that 
—t can be well approximated by the transverse momentum squared of the vector meson. 

In the present paper, the pQCD [J9Hl4[l and Regge [15 based approaches to vector-meson production 



are described. Then the measurements of elastic and proton-dissociative photoproduction of the 
p°, cj) and J/ifi mesons are presented. Finally, the results are compared with models and a summary 
of the conclusions is given. 



2 Models 

2.1 The pQCD approach — hard scale models 

In models based on pQCD 0-0], the diffractive photoproduction of a vector meson from a proton 
can be viewed in the proton rest frame as a three-step process: the photon fluctuates into a qq state; 



the qq pair scatters on the proton target; and the scattered qq pair becomes a vector meson. The 
qq fluctuation is described in terms of the photon wave- function derived from QCD. The interaction 
of the qq pair with the proton is mediated in the lowest order by the exchange of two gluons in 
a color-singlet state. In the leading logarithmic approximation (LLA), the process can also be 
described by the exchange of a gluon ladder |10|-p^|. To be calculable in pQCD, the interaction has 
to involve a hard scale or, in other words, has to involve small transverse distances. The transition 
of a qq pair into a meson is, however, a non-perturbative phenomenon and can only be described 
in terms of the meson wave-function derived from lattice calculations and sum rules 111 



The expected signatures of the perturbative regime in diffractive meson production are: 

• a fast rise of the diffractive cross sections with W, the available center-of-mass energy, due 
to the fast increase with decreasing x (Bjorken scaling variable) of the gluon density in the 
proton f9,10|; 



no variation with W of the ^-dependence of the cross section, i.e. no shrinkage of the 



diffractive peak 11]; 



approximate restoration of flavor-independent production, which is expected when the photon 
couples directly to the constituent quarks in the meson Q; 

production of light vector mesons in a helicity-zero state, independent of the initial photon 
helicity fl2|,[l4|], in the asymptotic limit of very large —t (W 2 S> — t 3> Aq CD ) . 



2.2 The Regge approach — modeling the soft interactions 

Regge phenomenology [15| has been successful in describing soft hadron-hadron interactions. In 



this approach, the interactions are described in terms of t-channel exchanges of Regge trajectories. 
In particular, diffractive processes are assumed to proceed through the exchange of the Pomeron 
trajectory. 

2.2.1 Regge factorization 



Regge factorization [15| is the assumption that Regge pole residues factorize into a contribution 
from each vertex. In other words, for diffractive vector-meson photoproduction, the properties of 
the interaction at the Pomeron-proton vertex should not depend on the properties of the Pomeron- 
vector-meson vertex. This hypothesis implies that the ratio of elastic to proton-dissociative vector- 
meson photoproduction, ^j(7P — ► Vp)/^( r yp - ► VN), should be the same for the three vector 
mesons under study in this paper. In the framework of VDM, these ratios should have the same 
values as in hadron-proton reactions. 



2.2.2 The Pomeron trajectory 

In general, the differential cross section for a two-body hadronic process, du/dt, can be expressed 
at high energies as 

^ = F{t)(W 2 f a ^- 2 \ (3) 



where F(t) is a function of t only and ot P (t) is the Pomeron trajectory. At lower W values, the 
exchange of a Reggeon trajectory should also be taken into account. 

By studying the W dependence of da/dt at fixed t, a p {t) can be determined directly. Usually, the 
trajectory is assumed to be linear, a p (t) = a p (0) + a' p t, but its form is not predicted by Regge 
theory. The early determinations of the Pomeron trajectory according to this procedure used data 
from pp elastic scattering []i~6|,|l7fl. 



Under the assumption that da/dt decreases exponentially (i.e. da/dt ~ exp (&£)), a' can also be 
determined from a study of the energy behavior of the exponential slope b. This method, however, 
is less direct and also depends on the t range over which the exponent is fitted. Determinations 



of a'jp based on this procedure [18, T^] from pp, Kp and irp elastic scattering yielded values in the 
range of 0.14-0.28 GeV -2 . However, the high precision pp ISR data at small —t showed that a! 
has a value 0.25 GeV -2 with a small uncertainty pQ] . The same analysis gave 0^(0) = 1.08. 

Studies of the elastic photoproduction of p° and 4> mesons have shown jj]] that these processes can 
be well described by the Regge phenomenology developed for soft hadron-hadron collisions. The 
steep energy behavior of the elastic J/ip photoproduction cross section at HERA, however, cannot 
be described in the Regge picture by a Pomeron trajectory with an intercept of 1.08 but requires a 
larger intercept. In addition, a direct determination of the Pomeron trajectory in a way similar to 
that described above has shown that the slope a'^ from elastic photoproduction of J/i/j is smaller 
than 0.25 GeV -2 [pi}] . These observations suggest that the Pomeron trajectory is not universal 
when a large scale, like a large mass, is involved. It is thus of interest to see whether the universality 
notion can be kept in soft interactions. Earlier attempts to determine the Pomeron trajectory from 



the elastic photoproduction of <fi mesons [22, 23 [ were not precise enough owing to the small span in 
the energy available. The present measurements at HERA, together with the existing lower-energy 
data, enable a more precise determination of the Pomeron trajectory and thus make it possible to 
test its universality. In this paper the Pomeron trajectory will be determined from the p° and 4> 
vector-meson data. The J/ip measurement of this analysis does not add significant information to 
the analysis done in pit] and therefore will not be considered here. 



3 Experiment 

The data used in the present analysis were collected with the ZEUS detector at HERA in 1995, 
when HERA collided positrons of energy E e = 27.5 GeV with protons of energy E p = 820 GeV. 
The data sample used in this analysis corresponds to an integrated luminosity of 1.98 pb _1 . 

A detailed description of the ZEUS detector can be found elsewhere [p3|j . A brief outline of the 
components which are most relevant for this analysis is given below. 

Charged particles are tracked by the central tracking detector (CTD), which operates in a magnetic 
field of 1.43 T provided by a thin superconducting coil. The CTD consists of 72 cylindrical drift 
chamber layers, organized in 9 superlayers covering the polar angle[|] region 15° < 6 < 164°. The 
transverse momentum resolution for full-length tracks is a{pt)/pt = 0.0058pi © 0.0065 © 0.0014/p 4 , 
with. p t in GeV ||. 



lr The ZEUS coordinates form a right-handed system with positive-.?? in the proton beam direction and a horizontal 
X-axis pointing towards the center of HERA. The nominal interaction point is at X — Y = Z — 0. The polar angle 
6 is defined with respect to the Z direction. 



The high resolution uranium-scintillator calorimeter (CAL) [^] consists of three parts: the forward 
(FCAL), the rear (RCAL) and the barrel (BCAL) calorimeters. Each part is subdivided transversely 
into towers and longitudinally into one electromagnetic section (EMC) and either one (in RCAL) or 
two (in BCAL and FCAL) hadronic sections (HAC). The smallest subdivision of the calorimeter is 
called a cell. The CAL energy resolutions, as measured under test beam conditions, are a(E)/E = 
0.18/VE for electrons and a(E)/E = 0.3b /VE for hadrons (E in GeV). 

The proton-remnant tagger (PRT1) is used to tag events in which the proton dissociates. It consists 
of two layers of scintillation counters perpendicular to the beam at Z = 5.15 m. The two layers are 
separated by a 2 mm thick lead absorber. Each layer is split into two halves along the Y-axis and 
each half is independently read out by a photomultiplier tube. The counters have an active area of 
30 x 26 cm 2 with a hole of 6.0 x 4.5 cm 2 at the center to accommodate the HERA beampipe. The 
pseudorapidity range covered by the PRT1 is 4.3 < rj < 5.8. 

The photoproduction tagger (PT) is a small electromagnetic calorimeter located at Z = —44 m, 
sensitive to 22-26 GeV positrons scattered under very small angles (less than a few mrad). The 
HERA positron beampipe has a 14 mm deep and 60 cm long indentation on the side facing the ring 
center (Fig. ||). The calorimeter is installed behind a 1 mm thick copper window in the beampipe. 
During beam injection and acceleration, a movable 10 cm thick lead shield is inserted in front of the 
PT. The detector consists of twelve 70 x 90 x 7 mm 3 tungsten plates interleaved with 3 mm thick 
scintillator layers. The light from the scintillator is read out from the bottom by a wavelength- 
shifter plate coupled through a plastic light-guide to a photomultiplier tube with a quartz window 
(Philips XP1911). The detector sensitive edge is about 28 mm from the positron beam. 

Additional scintillator strips are installed, after each of the first three tungsten plates, at depths 
corresponding to 2, 4 and 6 radiation lengths. These 8 mm wide vertical strips are connected to 
plastic light-guides coupled to three miniature photomultipliers (Hamamatsu R5600). Signals from 
the strips can be used to apply fiducial cuts and select well-contained electromagnetic showers. 

Test beam measurements demonstrated that for 1-5 GeV electrons hitting the calorimeter centrally, 
the energy resolution is a(E)/E = 0.25 j \J E '(GeV), and the calorimeter linearity is better than 1%. 
The energy measurement is used only at the trigger level for tagging photoproduction events. The 
tagger issues a trigger for events with an energy deposition above approximately 1 GeV. The low 
threshold ensures that the tagging efficiency is determined mainly by the geometric acceptance. 

e-beam 



Pb 

A 




PMT 
W/Sci calorimeter 



\\ (shield) 

Figure 1: Sketch of the layout of the photoproduction tagging detector. 

The luminosity is determined from the rate of the Bethe-Heitler bremsstrahlung process, ep — > ejp, 
where the photon is measured with a calorimeter (LUMI) located in the HERA tunnel downstream 



of the interaction point in the direction of the outgoing positron 27]. The acceptance of the LUMI 
calorimeter for Bethe-Heitler events is greater than 98%. 



4 Kinematics and cross sections 

The kinematics of the inclusive scattering of unpolarized positrons and protons are described by 
the squared positron-proton center-of-mass energy, s, and any two of the following variables: 

• Q 2 , the negative square of the exchanged photon's four-momentum; 

• y, the fraction of the positron energy transferred to the hadronic final state in the rest frame 
of the initial-state proton; 

• W 2 = ys + M 2 (l — y) — Q 2 , the squared center-of-mass energy of the photon-proton system 
(where M p is the proton mass); W 2 ~ ys in the case of photoproduction. 

For the exclusive reaction ep — ► eVp (V — > two charged particles) and the proton-dissociative 
process ep — ► eVN, t and the following additional variables are used (see Fig. g): 

• the angle, $, between the V production plane (which contains the momentum vectors of the 
virtual photon and the vector meson) and the positron scattering plane; 

• the polar and azimuthal angles, Oh and ip h , of the positively-charged decay particle in the V 
helicity frame. The polar angle, Oh, is defined as the angle between the direction opposite to 
that of the outgoing proton and the direction of the positively-charged decay particle. The 
azimuthal angle, </%, is the angle between the decay plane and the V production plane; 

• the mass, Mjy, of the diffr actively-produced state N in the proton-dissociative reaction. In 
the present analysis, however, it was not possible to measure this quantity directly and the 
Mn range covered was obtained from Monte Carlo simulations (see Sect. pi). 



'<Ph< 
V production plane / ) 



positron scattering plane 



V decay plane 




Figure 2: Illustration of the angles used to analyze the helicity states of the vector meson (for a 
decay into two particles, V — ► q + q~). 

Only the three-momenta of the decay particles were measured. Neither the momentum of the scat- 
tered positron, nor the $ angle, were measured. In such tagged photoproduction events, Q 2 ranges 
from the kinematic minimum, Q^ in ~ M 2 y 2 /(1 — y) ~ 10~ 9 GeV 2 , where M e is the electron mass, 
to a maximum value limited by the angular acceptance of the PT, (Jmax ~ ^E e E e i sin (0 max /2) ~ 



0.01 GeV 2 , where # max is the maximum scattering angle, and E e and E e / are the energies of the 
initial- and final-state positrons. The energy of the scattered positron is determined by the PT 
acceptance. Since the typical Q 2 is very small (the median Q 2 is approximately 7 x 10~ 6 GeV 2 ), 
it can be neglected in the reconstruction of the other kinematic variables. The photon-proton 
center-of-mass energy is given by 

W 2 ^2E P (E- PZ ), (4) 

where E p is the incoming proton energy, and E and pz are the energy and longitudinal momentum 
in the laboratory frame of the produced meson V, respectively; in this approximation, Mn, the 
meson mass, My, and its transverse momentum, pr, are assumed to be much smaller than W. The 
four-momentum transfer squared is given by 

t = (M 2 N - M 2 -Q 2 - M^f/AW 2 - (j>* - p*) 2 - 4p* lP * sm 2 (9*/2) = t - Ap^p* sin 2 (r /2), (5) 

where p*, p* are the magnitudes of the photon and meson momenta and 9* is the angle between 
them. The starred quantities are defined in the photon-proton center-of-mass system and to is the 
maximum t(9* = 0) value. For Q 2 <C — t,p^, My,M^ <C W 2 , t and t are given by 

t « t - p\ « -Ml{M 2 N - M 2 )/W 2 - p\. (6) 

Since the maximum —to value in the kinematic range covered by this analysis is small (—to ^ 
7 x 10~ 3 GeV 2 for V = p°, (f> and -t < 6 x 10~ 2 GeV 2 for V = J/ip) compared to the -t value, it 
can be neglected. The four-momentum transfer squared is then given by 

t » -p\ (7) 

The differential and integrated photoproduction cross sections for the processes jp — ► Vp{N) were 
obtained from the cross sections measured for the reactions ep — ► eVp(N). The cross sections are 
related by 



d 2 a ep 
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[i + (i-y) 2 


2(1 - y) Q 2 min 


dydQ 2 


2irQ 2 
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y Q 2 



a lp {W) 

ip(y,Q*) a lp {W), (8) 

where a is the fine structure constant and (p(y,Q 2 ) is the effective photon flux. A measured ep 
cross section can thus be transformed into a ^p cross section, 

a " p Jf^y,Q 2 )dydQ 2 * 7 ' W 

if a 7 p is independent of W (or y) in the studied region (where <& 7 is the integrated effective photon 
flux). 

5 Event selection 

Vector mesons were observed in the two-body decay channels p° — ► 7r + 7r _ , <fi — ► K + K~ , J/^ — v e + e _ . 
and J/^ — > /U + /U~ via the reconstruction of two oppositely-charged tracks in the CTD. The scat- 
tered positron was detected in the PT and the proton or its fragments escaped undetected or were 
tagged in either the FCAL or the PRT1. 

6 



5.1 Trigger 

ZEUS uses a three-level trigger system. At the first level, a coincidence between signals in the 
PT and a track candidate in the CTD was required. Additionally, it was required that the energy 
deposition in any of the FCAL towers closest to the beampipe should not exceed 1.25 GeV in order 
to suppress proton beam-related backgrounds. An upper limit of 1 GeV on the energy deposited in 
the LUMI was also imposed; this requirement suppressed events having a random coincidence with 
bremsstrahlung. The second and the third trigger levels were mainly used to reject non-diffractive 
backgrounds by requiring exactly two tracks pointing to the same vertex with a Z-coordinate 
compatible with that of the nominal position of the interaction point, \Vz\ < 60 cm. 

The trigger efficiency was studied using data samples selected by two independent triggers and was 
found to be about 90% for the elastic events used in this analysis and flat in all relevant kinematic 
variables. The efficiency for the proton-dissociative events was significantly lower, of the order of 
10%, due to the FCAL energy requirement. This requirement restricts the mass of the dissociative 
system to M N < 7 GeV. 



5.2 Offline requirements 

In the offline event selection the following conditions were imposed: 

• exactly two oppositely-charged tracks from a common vertex; 

• each track with pt > 0.15 GeV and \r/\ < 2.2; 



• the vertex coordinates in the range Vr = \ V^ + Vy < 0.7 cm and \Vz\ < 40 cm; 

• 85 < W < 105 GeV, thereby selecting a region of high and well understood tagging efficiency; 

• a maximum energy of 200 MeV (RCAL) and 250 MeV (BCAL) deposited in any calorimeter 
cell, with the exception of those matched to tracks; 

• energy deposition in any of the 8 FCAL towers closest to the beampipe less than 1.20 GeV. 

The number of events thus selected was 25446. 

The final identification of the vector-meson candidates was performed using cuts on the invariant 
mass of the track pairs measured in the CTD. The pion, kaon and electron masses were in turn 
assigned to the tracks, leading to the selection of 22823 p° (0.55 < M W7r < 1.20 GeV,MKK > 
1.075 GeV), 366 <p (0.99 < M KK < 1.06 GeV), and 120 J/t/j (2.98 < M ee < 3.13 GeV) meson 
candidates. The invariant mass resolution varies from 5 to 60 MeV depending on the meson type 
and the values of t and W. 

6 Acceptance corrections 



6.1 Monte Carlo generators 



The exclusive reaction ep — ► eVp was modeled using the DIPSI [28 1 Monte Carlo (MC) program. 
For the simulation of the reaction ep — > eVN, the EPSOFT Monte Carlo generator [p£j] was used. 
The Mn distribution at fixed t was reweighted to the dependence 

^^F « ^TTv (10) 

dM N 2 M P N {t) 

For p° and (f) production, the function (3{t) = i.i2e a64+0 - 3t2+0 - 04 * 3 + i.08e a854+0 - m2 was used; this 



function was found by parameterizing the ISR |30[ data on -^ w in pp single diffraction at large — t 



dM 



JV 



and Af^r < 0.1 W 2 plj . For J ftp production the average between the above f3(t) and (3(t) = 2.35 
(expected in Regge phenomenology for constant Oijp{t) = 1.175 J32J) was used. 

The t (and M^ for V=p°) distribution was reweighted in both generators so as to reproduce the 
measured distribution after reconstruction; the polar and azimuthal angular distributions of the 
decay particles in the helicity frame were also reweighted. 



6.2 Photoproduction tagging acceptance 

The geometric acceptance of the PT was simulated by a program which uses the HERA beam- 
transport matrices to track the positron through the HERA beamline. This program was tuned 
so as to reproduce the measured tagging efficiency, Ab, using Bethe- Heftier events, ep — > e^p. 
Ab was defined as the fraction of bremsstrahlung events with a photon measured in the LUMI 
photon detector when the PT fired the trigger. In Fig. ||a, Ab is shown as a function of the 
measured photon energy, £L, in the LUMI; the MC prediction is in reasonable agreement with the 
data. The acceptance for photoproduction events was determined using the geometric acceptance 
of the tagger and events generated according to the equivalent-photon approximation for positron 
scattering angles up to 3 mrad (for larger angles the acceptance is negligible); for these events, the 
positron was tracked through the HERA beamline. The photoproduction tagging efficiency, A t , 
was calculated as a function of the positron energy, E e i, (see Fig. |3]b). For the kinematic range used 
in this analysis, the average (cross-section weighted) PT acceptance was 70%. The systematic error 
was evaluated by changing the data sets used to tune the MC and by varying the photon energy 
scale and the position of the positron exit window in MC within their systematic uncertainties. 



6.3 Overall acceptance 

The generated events were processed through the same chain of selection and reconstruction pro- 
cedures as the data, thereby accounting for trigger as well as detector efficiencies (except for that 
of the PT) and smearing effects in the ZEUS detector. The reconstructed Monte Carlo events were 
then weighted with the function At(E e i) in order to account for the PT acceptance and efficiency. 

All measured distributions are well described by the Monte Carlo simulations. Some examples are 
displayed in Fig. [|. The overall acceptance in a given bin was then determined as the ratio of the 
number of accepted Monte Carlo events (weighted by At) to the number generated in the selected 
kinematic range. The acceptance, calculated in this manner, accounts for the geometric acceptance, 
the detector and reconstruction efficiencies, the detector resolution and the trigger efficiency. 
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Figure 3: a) The measured bremsstrahlung tagging efficiency (dots), Ab, as a function of the 
measured photon energy, E^, compared to the MC expectation (solid line), (b) Photoproduction 
tagging efficiency At as a function of the energy of the scattered positron E' e ; the shaded band 
represents the systematic uncertainty. 

Figure || shows the overall acceptance for elastic events as a function of t, ph and cos Oh- Inefficiencies 
at small —t, as well as the strong variation of the acceptance with ph for p° and cj> mesons, are 
mainly due to the relatively small opening angle between the decay particles, resulting in many 
very backward tracks (as can be deduced from the pseudorapidity distributions of tracks in Fig. H) 
outside the geometric acceptance of the CTD. 

The average acceptances are of the order of 30% for elastic and 5% for proton dissociative events, 
respectively. The much lower acceptance for the proton dissociative events is mainly due to the 
FCAL energy requirement imposed in the trigger. 



7 Backgrounds 



The dominant background sources are non-resonant 7r + 7r~ production for the p° analysis, p° pro- 
duction for the <j> analysis, and Bethe-Heitler 77 — ► e + e~/^ + /i~ production for the J /ip case. These 
backgrounds were statistically subtracted using the fits to invariant-mass distributions as described 



in Sect. 1C. 



The background due to inclusive photon diffractive dissociation, 7^ — » Xp, was studied using Monte 
Carlo simulations and minimum-bias data samples. It was found to be about 1-2% and was also 
subtracted in the fitting procedure. 

In the p° analysis, the backgrounds due to decays of lo and (j> mesons were found from MC studies 
to be negligible. The background due to non-diffractive events was also small and was neglected. 

On average, 10% of the diffractive events were rejected by the trigger due to accidental coincidences 
with bremsstrahlung events in the LUMI calorimeter. A correction was applied to account for this 



effect. A small fraction of events, below 1%, was selected due to fake photoproduction tagging 
from bremsstrahlung overlays (when the meson decay particles were measured by the CTD and the 
scattered positron was undetected, while the PT was hit instead by a bremsstrahlung positron); 
this effect was neglected. 



8 Separation of elastic and proton-dissociative processes 

The selected samples of diffractively-produced vector mesons are mixtures of elastic and proton- 
dissociative events. These two processes have been separated as a function of t on a statistical basis 
using the MC simulation. 

Proton-dissociative events were tagged by requiring a signal in one of the PRT1 counters above a 
threshold corresponding to the signal of a minimum ionizing particle. Alternatively, for systematic 
checks, energy depositions in the FCAL towers close to the beampipe were also used for tagging 
the proton-dissociative events. 

Since the non-diffractive backgrounds were negligible, the following relation was assumed: 

lyl pd,data _ Is pd,MC 
N p d,data N p d,MC 

where N p d,data and N p ^mc are the numbers of all accepted (but not necessarily tagged) proton- 
dissociative events, whereas N^ data and Np d MC are the numbers of tagged proton-dissociative 
events, in the data and proton-dissociative EPSOFT samples, respectively. Therefore, the fraction, 
Cd j of proton-dissociative events in the data was calculated from 

fi _ N p d,data _ ™pd,data ,^pd,MC 
Ndata Ndata Npd,MC 

where Ndata is the number of all observed events (elastic and proton-dissociative) in a given t bin 
in the data. In Fig. |6] the observed fraction of PRT1 tags, Rjj, in the data and in the proton- 
dissociative MC, as well as the estimated fraction of proton-dissociative events, Cjj, are displayed 
separately for the p°, 4> and J/tp samples. For the low —t region, the elastic process contributes 
a large fraction of the selected diffractive events. For —t > 1.0 GeV 2 , the contribution from the 
proton-dissociative process exceeds that from the elastic channel. 
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Figure 4: Comparison between the data and MC distributions of r] and pt of the positively charged 
track, and W, cos 6h, <Ph an d the meson pseudorapidity r)y. The three columns refer to the sum of 
the elastic and proton-dissociative p , and J/-0 samples, respectively. 
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Figure 5: Overall acceptance for elastically produced p°, <fr and J/ip mesons as a function of —t 
and the helicity variables <f>h and cos 6h- 
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Figure 6: Observed fraction of PRT1 tags, Rd, in the data (triangles with statistical error bars) 
and in the proton-dissociative MC (histogram), and the estimated fraction of proton-dissociative 
events in data, Cd, for the p° (a,b), (f> (c,d) and J/ij) (e,f) samples as a function of —t. The inner 
error bars indicate the statistical errors, the outer bars the statistical and systematic uncertainties 
added in quadrature. The shaded bands in b), d) and f) represent the size of the correlated errors 
due to modeling of the proton dissociation in the Monte Carlo. 



13 



9 Systematic uncertainties 

The systematic uncertainties were subdivided into those related to the selection procedure and the 
detector simulation, and those reflecting the uncertainty in the model used for the Monte Carlo 
generator. 

The effects of the following changes in the selection procedure were checked: 

1) the FCAL (instead of PRT1) was used to tag the proton-dissociative events; 

2) the matching procedure between the tracks and energy depositions in CAL was varied by 
changing the matching distance between the track and CAL cluster and the energy require- 
ment used to define the CAL cluster; 

3) the minimum value of the track transverse momentum was varied between 100 and 200 MeV; 

4) the limit on the track \rj\ was varied between 2.1 and 2.3; 

5) the requirements on the vertex position were varied by amount corresponding approximately 
to the resolution, i.e. ± 5 cm in Z and ± 0.25 cm in the radial direction; 

6) the selected invariant mass region was changed by widening and narrowing it by amounts 
corresponding to approximately the appropriate mass resolution; 

7) the threshold values of all PRT1 counters were increased by 100%. 

The first check resulted in 5-20% changes of the measured cross sections. Checks 2-5 resulted in 
3-5% changes, and the last two checks had negligible effect. 

In the MC simulation, the positions of some beamline elements (eg. the position of the synchrotron 
collimator jaws) were varied within their uncertainties, resulting in a 5-15% change in the cross 
section. 

In order to estimate the systematic uncertainty from the uncertainties in the parameters assumed 
in the EPSOFT generator, the following modifications were made: 



• the shape of (5{t) was varied within the uncertainties of its measurement [31|; for p° and cf) 
production, this corresponds to (3{t) = 2.24e a6<+a3t2+0 - 04 ' 3 and (3(t) = 2.16e a85t+0 - 11<2 ; for 
J/V> production, to f3(t) = i.i2e a6 '+ a3 ' 2+a04t3 + i.08e a85t+a11 ' 2 and (5{t) = 2.35. These 
variations significantly changed the acceptance corrections, by up to 10% at low —t and up 
to 70% for the proton-dissociative sample at the highest —t; 

• multiplicity distributions of the decay particles of the dissociative system N were varied within 



the uncertainties of their measurement [33|. This resulted in cross section variations of < 5%. 



Additionally, the re-weighting of other MC distributions (of the decay particle angles or invariant 
masses, for example) was performed in the range allowed by maintaining satisfactory agreement 
between data and Monte Carlo. The effect on the cross sections was < 3%. In the determination of 
the p° spin-density matrix elements, the difference between the nominal method (see Sect. 10.1.3 ) 
and the method of moments was taken as an additional uncertainty. 

The overall normalization error due to the photoproduction tagging uncertainty was ±15%. 
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10 Results 



The differential cross sections da/dt for elastic and proton-dissociative photoproduction of V were 
evaluated in each bin of t as: 



da 7P ^ V p _N-(l-C D )- C r 



dt 



A- €■$-,■ At- C br 



(11) 



da. 



7P— >VN 

dt 



N-C D -C r 



A-C-^^-At-Cbr' 



(12) 



where iV is the number of observed vector-meson candidates in bin At after all selection cuts, Cd is 
the estimated fraction of the proton-dissociative events in the bin, C res is the resonant contribution 
in the bin, C& r is the branching ratio of the vector-meson decay mode considered, A is the overall 
acceptance in the bin, C the integrated luminosity, and <3? 7 is the integrated effective photon flux 
(see Eq. |9|). In the kinematic region Q 2 < (E e 6 max ) 2 (l — y), where max = 3 mrad (see Sect. 6.2), 
and 85 < W < 105 GeV, $ 7 = 0.0121. The branching ratios for the p° -> vr+vr-, <p -> K+K~, 
J/tp — ► e + e~ or J/tp — » fi + fj,~ decay modes were taken as 1, 0.5 and 0.12, respectively. Effects due 



to QED radiation, estimated to be smaller than 2% [34|, were neglected. 



A similar procedure was used to evaluate the differential cross section da/dM^- 

For the proton-dissociative reaction, cross sections are extrapolated to M^ = O.liy 2 using the 



EPSOFT MC, as modified by Eq. 1C. 



10.1 p° photoproduction 



10.1.1 Resonance mass shape 



The differential cross sections da jdM^^ were fitted using a parameterization inspired by the Soding 
model p5[, where the p-wave relativistic Breit-Wigner (BW) shape is distorted by the interference 
with non-resonant tttt production: 



da 

~d~MZ., 



A 1 



y/M, 



T7T lVIpL p 



ML - M 2 + iM D T 



+ B/A 



pj- p 



+ f 



(13) 



where M„ is the p° mass, T p is the momentum-dependent width 



I\,(M fl 



go 



Mp 

M^ 



(14) 



To is the width of the p°, q is the 7r momentum in the tttt rest frame and go is the value of q at 
M nn = M p . The non-resonant amplitude (taken to be M^-independent) is denoted by B, and A is 
the normalization factor of the resonant amplitude. Additionally, another term, /, was introduced 
to account for the background from reactions with photon diffractive dissociation. The term / was 
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assumed to be linear in M nir , / oc (1 + 1.5M„- 



with M n7T in GeV. Alternatively the following 



parameterization, proposed by Ross and Stodolsky fl36l, was used: 



da 

dMir- 



A 



M« 



;M p T p 



(M; 



7T7T 



Mty + MX 



Mp 



+ f 



(15) 



where T p is given by Eq. 14 and / has the same form as described above. 



The x 2 /NDF for all the fits is satisfactory. The fitted values of M p , To and / do not depend on 
the prescription used to parameterize the mass distribution. The mass, M p , and width, Tq, for 
all t values are compatible with the values of the Particle Data Group (PDG) |37|| . The level of 
background under the p° peak, as given by the integral of the function /, is about 1-2%. The mass 
distributions for elastic and proton-dissociative p° production are shown in Fig. together with 
the results of the fits using Eq. ^. The mass resolution varies between 20 and 60 MeV for Af ff7r 
between 0.5 and 1.2 GeV, respectively. 
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Figure 7: The differential cross sections da/dM n7T for several t ranges. The points represent the 



data and the curves indicate the result of the fits with Eq. 13. The dashed curves represent the 
resonant contribution, the dotted curves the non-resonant contribution and the dot-dashed curves 
the contribution from the interference term. The solid curves are the sum. Only the statistical 
errors are shown. The circles and squares correspond to the elastic and the proton-dissociative 
samples, respectively. 



16 



The fits to the mass distributions were repeated with M p and Tq fixed to the PDG values. The 
results for B/A and n are shown as a function of —t in Fig. B; they match with our earlier measure- 
ment for the low —t region at (W) = 70 GeV ||. Both B/A and n decrease with —t, indicating 
that the resonance shape distortion decreases. This decrease is much faster with —t than that with 



(16) 



The dependence of the ratio B/A on t was parameterized as 

B/A = ke bst . 



The fit to the data yields k = -0.86 ± 0.52 (stat.) GeV" 1 / 2 and b s = 1.6 ± 0.9 (stat.) GeV" 2 
for the elastic reaction, and k = -0.71 ± 0.21 (stat.) GeV~ 1/2 and b s = 0.9 ± 0.4 (stat.) GeV -2 
for the proton-dissociative process. The similar magnitudes and t dependences of the B/A ratio 
for the elastic and the proton-dissociative processes indicate that the data are consistent with the 
hypothesis of factorization of the diffractive vertices. 

The p° production cross section was evaluated in the mass range 2M n < M n7T < M p + 5IV The 
extrapolation from the measured range 0.55 < M wn < 1.2 GeV was made using the results of the 
fit to Eq. O. 
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Figure 8: The ratio B/A from Eq. 13 and the parameter n from Eq. 15, as a function of —t. 
The inner error bars indicate the statistical errors, the outer ones the statistical and systematic 
uncertainties added in quadrature. The solid circles and squares correspond to the elastic and 
the proton-dissociative samples, respectively. The open circles correspond to the low —t ZEUS 
results ||. The solid curves are fits of the form ke st (see text). 



17 



10.1.2 t— distributions 

The differential cross sections, da/dt, for elastic and proton-dissociative (for Mfj < 0.1W 2 ) p° 
photoproduction are plotted in Fig. ||a. Both samples exhibit an exponential drop with increas- 
ing —t, with the cross section for the elastic process falling off more steeply than that for the 
proton-dissociative process. A fit with the function ^4exp(6t) in the range 0.4 < — t < 1.2 GeV 2 
for the elastic process gives b e i = 6.0 ±0.3 (stat.) j^J (sys.) ± 0.4 (mod.) GeV -2 and b p d = 2.4 
±0.2 (stat.) 1q' 2 (sys.) ±0.3 (mod.) GeV~ 2 for the proton-dissociative sample, where (mod.) rep- 
resents the uncertainty due to the modeling of the proton dissociation in the EPSOFT Monte 
Carlo. 

The ratio of the elastic to the proton-dissociative cross sections, dcr(7p — ► /?°p)/dcr(7p — > p°N), in 
a given t interval is shown in Fig. |9p. This ratio falls off rapidly with —t from about unity at —t ~ 
0.4 GeV 2 to about 10~ 2 for — t > 1 GeV 2 . By comparison, this ratio has the value 2.0 ± 0.2 (stat.) 
± 0.7 (syst.) for —t < 0.5 GeV 2 at W = 70 GeV ||. An exponential fit to the ratio gives a value 
for the slope Ab = 3.4 ± 0.3 (stat.) _ ' 3 (sys.) ±0.2 (mod.) GeV -2 , which is similar to that found 
previously M. 



Figure 10 shows the differential cross section for the elastic reaction together with results obtained 



in our earlier studies of the low —t region j|, 39 1. The data at low — t were measured at a somewhat 
lower W ((W) = 70 GeV) and have been rescaled to the average W of the present analysis ((W) = 



94 GeV), assuming the Pomeron trajectory as measured in this analysis (see Sect. 11.4). The large 
—t data match well with the low —t results. 



The data in Fig. 10 cannot be described by the expression Aex.p(bt) over the whole t range. The 
slope b e i (of about 10 GeV ||) in the low —t region is larger than that in the large —t region, in 
contrast to the case for Ab, which is the same within errors in both t regions, as well as at higher 

iol. 



10.1.3 Decay angular distributions 

The angular distributions of the decay pions were used to determine some of the p° spin-density 
matrix elements using the invariant mass selection 0.55 < M irK < 1.2 GeV. The direction of the 
virtual photon was approximated by that of the incoming positron. In this measurement, the 
three-dimensional angular distribution has been averaged over the azimuthal angle between the 
positron scattering plane and the p° production plane, $, and thus no longer distinguishes the 
photon helicity states ± 1. The normalized two-dimensional decay angular distribution can be 
written as Q 

W(cos9 h , fa) = £{{h(l- r°ot) + i(3r ° 4 - 1) cos 2 9 h 

-V2Re[r%] sin 26 h cos <j) h -r^ sin 2 9 h cos 2<f> h } , (17) 

where the spin-density matrix element Tqq represents the probability that the produced p° has 
helicity 0; the element r^g is related to the interference between the helicity non-flip and single-flip 
amplitudes; rji 1 is related to the interference between the non-flip and double-flip amplitudes. If 

r.04 ^04 OT1 j T> Q U)4l 



s-channel helicity conservation (SCHC) E^] holds, Tqq, r^_ 1 and Refr^] should be zero. 

The parameters Tqq, r^ 4 _ 1 and Re[r° 4 ] were obtained by minimizing the difference between the two- 
dimensional (cosOh, <frh) angular distributions of the data and those of the simulated events, which 
were re- weighted according to Eq. 17. A binned x 2_me thod was used. 
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Figure 9: a) The differential cross sections du/dtfor elastic (circles) and proton-dissociative 
(squares) p° photoproduction. The solid lines represent the results of the fit with the function 
Ae . The normalization error of 15% is not shown, b) The ratio of the elastic to the proton- 
dissociative cross sections shown in a). The inner error bars indicate the statistical errors, the 
outer bars the statistical and systematic uncertainties added in quadrature. The shaded bands 
represent the size of the correlated errors due to the modeling of the proton dissociation in the 
Monte Carlo. 



The three spin-density matrix elements are shown in Fig. |11 
dissociative p° production as a function of 



separately for elastic and proton- 
tends to be 



-t. In the t range of this analysis, Refr^Q 
non-zero and positive, while rj_i tends to be non-zero and negative. This is an indication for small 
deviations from SCHC giving rise to non-zero helicity single-flip and double-flip amplitudes. The 
element Tqq is measured to be zero within errors. 

The present results are not corrected for the non-resonant ir + ir~ production. Hence, strictly speak- 
ing, they only apply to the reaction jp(N) — > Tr + ir~p(N) within the quoted M nn range. However, 
in the previous section it was shown that the non-resonant contribution is small and decreases 
rapidly with —t. In order to assess the sensitivity of the data to changes in the selected M n7T range, 
the events with M n7T < M p and M n7T > M p were analyzed separately. No statistically significant 
effect was observed. Data at lower —t (< —t >w0.1 GeV 2 ) from ZEUS || show no evidence for a 
violation of SCHC. However, lower energy (< W >« 4.3 GeV) elastic photoproduction data |43|| , 
while showing no SCHC violation for —t < 0.2 GeV 2 , yield positive values for Refr^] and nega- 
tive values for rj ) l 1 for —t > 0.2 GeV 2 . Recent HERA measurements, both at low Q 2 (0.25-0.85 
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Figure 10: The differential cross sections, da/dt, for elastic p° photoproduction. The solid circles 
are those shown in Fig. pt the open circles and crosses are from earlier ZEUS data 0,^9|. The shaded 
band represents the size of the correlated errors due to the modeling of the proton dissociation in 
the Monte Carlo. The normalization error of 15% is not shown. 



GeV 2 ) Q and higher Q 2 (greater than 1 GeV 2 ) |44],|I|], also report similar values, at < — t >? 
0.1-0.2 GeV 2 , although only the low Q 2 results are significantly non-zero. 



10.2 photoproduction 



10.2.1 Mass distribution 



The invariant mass of the two charged decay products of the (j> is displayed in Fig. 12 for different 



—t ranges for two samples of data: for all events and for those that have a signal in PRT1 (proton 
dissociation). The invariant mass was computed assuming that the two charged particles were 
kaons. The lines are fits to a Breit-Wigner function, convoluted with a Gaussian resolution function, 
and a function describing the background. The background, due mainly to diffractive p° events for 
which the two pions are assigned the kaon mass, was assumed to have the form ~ {Mkk — 2-Mr-) , 
where A is a parameter determined by the fit. For the determination of the background contribution 
in each t region, the mass and width of the cj) were fixed to the PDG values. The fit was performed 
over a broad mass range 0.99 < Mkk < 1-13 GeV to give a better estimate of the background. 
A narrower mass range 0.99 < Mkk < 1-06 GeV was then used to select (j) candidates and the 
background contribution in this range was subtracted. The mass resolution is about 10 MeV, 
consistent with the Monte Carlo simulations. The background contribution decreases with —t, 
from about 23% at -t = 0.5 GeV 2 to 10% above -t = 1.2 GeV 2 . 
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Figure 11: The spin-density matrix elements ?qq, Re[r^] and r®_ 1 as a function of —t obtained by 



fitting Eq. |17] to the data. The squares correspond to the proton-dissociative sample and the solid 
circles to the elastic sample. The inner error bars indicate the statistical errors, the outer bars the 
statistical and systematic uncertainties added in quadrature. 



10.2.2 i-distributions 



The differential cross sections dcr/dt for the reactions jp 



and 7 p -► cj)N {M% < 0.1W 2 ) 



were determined following the separation procedure described in Sect. Q the results are shown in 
Fig. 13a. As in the p° case, both differential cross sections decrease exponentially with increasing 
—t. A fit with a function of the form dcr/dt = Aex.p(bt) in the range 0.4 < — t < 1.2 GeV 2 yields b e i 
= 6.3 ±0.7 (stat.)±0.6 (sys.)±0.3 (mod.) GeV" 2 and b pd = 2.1 ±0.5 (stat.)±0.3 (sys.)±0.4 (mod.) 
GeV" 2 . These values are the same within errors as those obtained for the p°. Note that while 



b p d > W el in the region -t < 0.3 GeV 2 |||, b p d ss W el in the region 0.4 < -t < 1.5 GeV 2 . There are 



no data at low —t to make a similar comparison in the case of the proton-dissociative process, but 
for the larger —t range, b p d as b d . This indicates that the mass difference between the p° and the 
<j) is not important in the large —t region. 



Fig. 13b shows the ratio of the cross sections for the elastic and the proton-dissociative processes 
in a given t interval as function of —t. The ratio decreases rapidly from unity at —t ~ 0.5 
GeV 2 to w 10" 2 at -t « 1.3 GeV 2 , as in the p° case. A fit to the ratio gives A6 = 3.0 ± 0.8 
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Figure 12: The Mkk distributions in four —t bins. The left-hand plots are for all events, while the 
right-hand plots are for PRTl-tagged events. The points represent the data and the solid curves 
indicate the result of the fit discussed in the text. The dashed curves represent the background 
contribution. Only the statistical errors are shown. 



(stat.) ^q'q (sys.) 1q*5 (mod.) GeV 2 . This value is the same within errors as that obtained for the 
p°, and is thus independent of the type of the vector meson produced at the photon vertex (see 
Sect. [TOI) . 

In Fig. [T| the results for da/dt for elastic <fi photoproduction are displayed, together with ZEUS 
measurements at small —t f46|| . The low —t data have been rescaled to the present W value by 
assuming the Pomeron trajectory as measured in this analysis (see Sect. 11. 4j ); both data sets are 
plotted with statistical errors only. The results are consistent. 



22 



ZEUS 1995 




0.25 



0.5 0.75 



■e- 
t ' 



1.25 1.5 1.75 2 2.25 

-t (GeV 2 ) 



Q. 



to 1 



a. 
-e- 

t 
a. 

$10 

b 
-O 




b) 



0.2 



0.4 



0.6 




1.2 1.4 

-t (GeV 2 ) 



Figure 13: a) The differential cross sections da/dt for elastic (circles) and proton-dissociative 
(squares) <j) photoproduction. The solid lines represent the results of the fit with the function Ae . 
The normalization error of 15% is not shown, b) The ratio of the elastic to proton-dissociative 
cross sections. The inner error bars indicate the statistical errors, the outer bars the statistical 
and systematic uncertainties added in quadrature. The shaded bands represent the size of the 
correlated errors due to the modeling of the proton dissociation in the Monte Carlo. 
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Figure 14: The differential cross sections da/dt for elastic (ft photoproduction. The solid circles 
correspond to the data of this analysis and the open circles are from an earlier ZEUS result |h| . 
The inner error bars indicate the statistical errors, the outer bars the statistical and systematic 
uncertainties added in quadrature. The shaded band represents the size of the correlated errors 
due to the modeling of the proton dissociation in the Monte Carlo. The normalization error of 15% 
is not shown. 
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10.3 J/ip photoproduction 



10.3.1 Mass distribution 



Figure |15| shows the two body invariant mass distribution in the region 2.4-3.6 GeV, for the whole 
sample in two —t ranges and for the PRTl-tagged sample in one —t bin. No lepton identification 
was performed for this sample. Although the events therefore represent a sum of e + e~ and /x + /U~ 
final states the electron mass was always assumed for each of the two charged particles^]- A narrow 
peak is observed around 3.1 GeV. Fits were made to the invariant mass distribution using the sum of 
a Gaussian (for the \i case) , a Gaussian (modified by the energy loss spectra due to bremsstrahlung 
in dead material before the CTD for the e case), and a background which linearly decreases with 
mass. Equal numbers of e pairs and \x pairs were assumed. The fitted J/ip mass is in good 
agreement with the PDG value. The background in the J/ip mass region (2.87-3.13 GeV) is (20 ± 
5)%, independent of t in the range under study {—t < 3 GeV 2 ). 
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Figure 15: The M ee distributions in two —t bins for the whole sample (a,b) and in one bin for the 
PRT1 tagged events (c). The points represents the data and the curves indicate the result of the 
fit discussed in the text. The dashed curve represents the background contribution. 



2 Assuming the muon mass does not change any of the results. 
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10.3.2 t— distributions 

The differential cross sections da/dt for the elastic and the proton-dissociative J/ip photoproduc- 
tion reactions are shown in Fig. |j~6|a. In contrast to the p° and (f) cases, the high mass of the 
J /if) and the resulting large opening angle of the decay particles result in significant acceptance 
in the low —t region. However, in the case of the proton-dissociative process, the effect of the 
minimum kinematically allowed —t, —to ~ 0.06 GeV 2 (see Eq. g), prevents the measurement of 
da/dt in the low —t bins. The elastic differential cross section falls exponentially, with bj = 
4.0 ±1.2(stat.)l 11 (sys.)l ' 6 (mod.) GeV -2 , while for the proton-dissociative reaction b d = 0.7 

±0.4(stat.)±0.2 (sys.)_g' 3 (mod.) GeV -2 . The slope Ab agrees, within errors, with the values for 
the p° and (j>. 



Fig. |l6| b shows the ratio of the elastic to the proton-dissociative differential cross sections. This 
ratio also falls from a value of about 1 at low —t to a value of about 0.1 for — t > 1 GeV 2 , similar 
to the p° and the 4> cases. 
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Figure 16: a) The differential cross sections da/dt for elastic (circles) and proton-dissociative 
(squares) J/ip photoproduction. The solid lines represent the results of the fit with the function 
Ae . The normalization error of 15% is not shown, b) The ratio of the elastic to the proton- 
dissociative cross sections. The inner error bars indicate the statistical errors, the outer bars the 
statistical and systematic uncertainties added in quadrature. The shaded bands represent correlated 
errors due to the modeling of the proton dissociation in the Monte Carlo. 
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The present measurements of da/dt for the reaction 7p — ► J/i^p are shown in Fig. [17], together 
with our earlier untagged (i.e. Q 2 ~ 0) photoproduction data |P| ((W) = 90 GeV). The data sets 
agree within errors. 
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Figure 17: The differential cross section da/dt for elastic J/ip photoproduction. The solid circles 
correspond to the present measurement, the open circles to the published ZEUS untagged photo- 
production results |]4q] . The shaded band represents correlated errors due to the modeling of the 
proton dissociation in the Monte Carlo. The normalization error of 15% is not shown. 



11 Comparisons of data to models 



11.1 Cross section comparisons with the pQCD-based models 



In order to compare the data to pQCD predictions for light |l | and heavy | ij mesons, the cross 
sections were redetermined in the region of validity of the model calculation, viz. Mjy < 0.01W 2 , 
using the EPSOFT MC simulation. The measured proton-dissociative cross sections for the three 
vector mesons are shown in Fig. [l^. 

The calculations [Q for the production of the p° and (j) mesons were performed at lowest order in 
as- Only helicity non-flip and single- flip amplitudes were taken into account. The non-perturbative 



effects were simulated by introducing effective quark masses. In Fig. 18a and 18b, the solid (dotted) 
curves represent results for a quark mass of 300 (200) MeV. For the calculation of the p° cross- 
section, the sub-asymptotic p° wave-function was used, whereas for the cj) production the asymptotic 
one [12, 14] was employed. The contribution of the perturbative cross section, represented by the 
dashed lines in Fig. |l8| a and 18b, is well below the p° and <p data. This observation, together with 
the helicity analysis (see Sect. 10.1.3; ) in which no significant production of p° mesons with helicity 
is observed, implies that these perturbative calculations are not applicable in this regime. 
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The situation is different in the J/tf) case. The perturbative QCD prediction in the LLA [|l 
expansion in terms of In (W 2 /Wq) (equivalent to In 1/x) compares satisfactorily with the data 



(solid line in Fig. |l8|c). A value of 015=0. 2 and an energy scale Wq = 1 GeV have been used in this 
calculation. It should be noted, however, that the uncertainties due to the choice of as (±10% - 
dotted lines) and the Wo scale (0.2 < Wq < 5 GeV 2 - dashed-dotted lines) are significant. 
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Figure 18: Comparison of the measured differential cross sections dcr/difor proton-dissociative 
vector-meson production for Mfj < 0.01VF 2 with the pQCD-based models for a) p° , b) 4> |0| and 
c) J/ip [ i^]. The inner error bars indicate the statistical errors, the outer bars the statistical and 
systematic uncertainties added in quadrature. The shaded bands represent correlated errors due 
to the modeling of the proton dissociation in the Monte Carlo. The normalization error of 15% is 
not shown. For the description of the curves see text. 



11.2 Ratios of cross sections for vector-meson photoproduction 

Flavor independence predicts that the ratio of the production cross sections of 4>:p° should be 
2:9 and that of J/ip:p° should be 8:9. These predictions are in striking disagreement with the 



28 



previously published low-Q 2 data: at W=70 GeV, the ratios 4>:p° = 0.065 ± 0.013 (47| and 
J/ip:p° = 0.00294 ± 0.00074 [03] were measured. In contrast, in the DIS kinematic region, for 
Q 2 « 12 GeV 2 , the ratios </>:p°=0.18±0.05 @ and 0.19±0.04 g| were obtained at W ss 100 GeV. 
This analysis gives the ratio J/V>:p°0.64±0.35 for Q 2 = 10 GeV 2 and 1.3±0.5 for Q 2 = 20 GeV 2 @ 
at W ~ 100 GeV. These results suggest that flavor independence may hold at large Q 2 . 

The ratios (j) '■ P° and J/ip '■ P° of the differential cross sections for elastic reactions from this analysis 



are plotted in Fig. |19|a and 19c as a function of —t. For both ratios, the point at t = was obtained 



by using the ratios of values of the total elastic cross section and rescaling them by the ratios of 
the corresponding slopes, b, of the differential cross sections, using the relation da/dt{t = 0) = bo. 
An increase of the eft : p° ratio up to — t ~ 1 GeV 2 is observed, approaching the expected 2:9 ratio. 
Although the J/ifi : p° ratio increases quickly up to — t ~ 0.5 GeV 2 , it remains more than an order 



of magnitude smaller than the 8:9 expectation. A non-perturbative QCD model 
describes the ratios for the elastic production at —t > 0.3 GeV 2 . 



successfully 



The ratios for the proton-dissociative reactions are shown in Fig. |19[ b and |19|d. There are no data 
at t = for this process. Both ratios are consistent with the corresponding elastic results. At 
—t ~ 2.2 GeV 2 , the J ftp : p° ratio is still well below the 8:9 expectation. 
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Figure 19: The ratios of the cross sections d<r/di for to p° (a and c) and J/ip to p° (b and d) 
for elastic and proton-dissociative photoproduction. The inner error bars indicate the statistical 
errors, the outer bars the statistical and systematic uncertainties added in quadrature. The ratios 
at t = were obtained from the data [47], 48] as explained in the text. The solid curves represent 
predictions of a non-perturbative QCD model |5lj for elastic photoproduction. The dashed lines 
correspond to the expectations of flavor independence. 
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11.3 Test of Regge factorization hypothesis 



The hypothesis of Regge factorization |52J implies that the ratio of the elastic to proton-dissociative 
differential cross sections should be independent of the type of vector meson produced at the 
photon vertex. Furthermore, this ratio should be the same as for hadron-proton elastic and proton- 
dissociative processes. 

Several results discussed in previous sections are consistent with this hypothesis: the decrease in the 
7T7T resonance-shape distortion with increasing —t is the same for the elastic and proton-dissociation 
channels; the p° spin-density matrix elements for both channels agree; and the difference (A6) 
between the slopes of the ^distributions for elastic and proton-dissociation processes is independent 
of the type of vector meson (/?°, cj>, J/ip) produced at the photon vertex. 

To test this hypothesis further, the ratios -SridP - > Vp)/-Tf{lP ~^ VN) for the three vector mesons 



V = p°,(j), J/tp are shown in Fig. 2C. As can be seen, the ratio agrees within errors for each of 
the three vector mesons. In the same figure, the ratio -§-{pp — ► PP)/-§;{pp -^ pN) at center of 
mass energies of y^ = 23.4 and 38.3 GeV Q| is also shown. The pp proton-dissociative reaction 
is defined for M^ < 0.05s. The ratios for the pp reactions are in agreement with those of the 
vector-meson photoproduction. These observations confirm the factorization hypothesis. 



In a naive additive quark model [53], the ratio shown in Fig. 20 measures the proton form factor. 
In such a model, the differential cross sections can be related to the form factor of the hadrons 
involved in the scattering process. In the framework of the VDM, this relation can also be used for 
the vector-meson photoproduction reactions discussed above. The elastic and proton-dissociative 
processes can be expressed as 



^ ^ F? F} \A\\ 

(*V* W> - (18) 

where A is the amplitude describing the constituent interaction, the form factor of the proton 
is indicated by F p and that of the vector meson by Fy- Thus, the ratio of the elastic and the 
proton-dissociative cross sections gives the proton form factor: 



da 
dt 



(7P -> Vp) _ _ 2 _ D 



where the electromagnetic dipole expression for the proton form factor has been used and D and 
to are free parameters. 

Equation ^ was fitted to both the photoproduction and pp data and the result is shown as the 
curve in Fig. |o|. The best fit yielded the parameters D = 7.7 ± 0.4 and m 2 = 0.47 ± 0.01 GeV 2 , 
and gives a fair description of the data. Note that the mass scale parameter obtained by this fit 
is smaller than the value of 0.71 GeV 2 measured from electron-proton elastic scattering.. Fitting 
only to the ZEUS photoproduction data yields D = 11.2 ± 8.0 and m 2 = 0.43 ± 0.13 GeV 2 . 

11.4 Determination of the Pomeron trajectory 

The measurement of the variation of the energy dependence of the elastic cross section with mo- 
mentum transfer t yields a direct determination of the Pomeron trajectory, as shown in Eq. ||. Such 
an analysis is presented for elastic p° and <fi photoproduction. 
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• V=p° ZEUS 7p 
o V=* ZEUS rp 
V=J/V ZEUS yp 

▼ p p Vs=23.4 GeV : 
' p p Vs=38.3 GeV 




-t (GeV 2 ) 

Figure 20: The ratio of the elastic to proton-dissociative differential cross sections as a function of 
—t for vector-meson photoproduction, together with data from pp reactions |3(| at y/s = 23.4 and 
38.3 GeV. The curve is the result of a combined fit to all the data as explained in the text. 



Since the reaction jp — > p°p is dominated by Pomeron exchange only at high energies, d<r/dt 
measurements at very low W cannot be used in this analysis. Therefore, the only fixed target 
experiment that can be used is the OMEGA experiment (WA4) [54|, which measured elastic pho- 
toproduction at W = 8.2 and 10.1 GeV in the t range 0.06 < —t < 1 GeV 2 . The measurement of 
the HI collaboration |55| at W=55 GeV, the earlier measurements of the ZEUS collaboration at 
W=71.2 GeV f§ and at W=73 GeV ^, and the present data at W=94 GeV are also used. 



The da/dt data used in the determination of a P (t) are presented in Fig. 21 in 12 t bins in the 
range < — t < 0.95 GeV 2 . The errors are the statistical and systematic uncertainties combined in 
quadrature. The line in each t bin is the result of a fit of the form W n , where n = 4a p (t) — 4. The 
resulting values of a F (t) are plotted in Fig. ^ as a function of t. A linear fit to the data yields 



a P (t) = (1.096 ± 0.021) + (0.125 ± 0.038)*, 



(20) 



and is plotted as a full line in the figure. The quality of these fits is acceptable. The dashed line 
is the Pomeron trajectory 1.0808 + 0.25 *, as determined by Donnachie and Landshoff (DL) |5T 
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The resulting intercept, 0:^(0), is in excellent agreement with that of DL; however, the slope, oi , 
is smaller in the present determination. 
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Figure 21: Cross sections for exclusive p° production from ZEUS, HI and OMEGA |j, 39, 54, 55] at 
fixed — t values as a function of W. The error bars show the statistical and systematic uncertainties 
added in quadrature. The lines correspond to the results of the fits to da /At oc (W 2 ) 2a jp^'^ 2 . 
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Figure 22: Determination of the Pomeron trajectory from the reaction 7p — » p°p. The dots are 



the values of the trajectory at a given — £ as determined from Fig. 21 and the solid line is the 



result of a linear fit to these values. The Pomeron trajectory as determined by DL [En] is shown 
for comparison as a dashed line. 
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The elastic photoproduction of <f> mesons is a good reaction with which to study the properties of 



the Pomeron, since this is the only trajectory that can be exchanged [57 1, assuming the ^ to be a 
pure ss state. This allows the use of data at very low W. 



The da/dt data used for the trajectory determination include the following: iy=2.64-3.60 GeV |5 
W=2.8jK$, ^=2.81-4.28 GeV g§, VF=3.59-4.21 GeV §§, VF=4.73-5.85 GeV @, VF=12.89 
GeV pl| , W=70 GeV [fl(| and the present measurement at VF=94 GeV. These data points are 
displayed in Fig. 23 for 11 t values in the range < — t < 1.4 GeV 2 . The lines are the results of 
fits of the function W n to the data. 
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Figure 23: Cross sections for exclusive <f> production from ZEUS and low energy measurements 



43 , 46, 58 31 at fixed —t values as a function of W . The error bars show the statistical and 
systematic uncertainties added in quadrature. The lines correspond to the results of the fits to 

da/dt (x{W 2 ) 2a ip^- 2 . 
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The resulting values of the trajectory a p {t) are shown in Fig. 24 as a function of t. Assuming a 
linear trajectory, 



a P (t) = (1.081 ± 0.010) + (0.158 ± 0.028)* 



(21) 



is obtained. Again the quality of the fits is generally acceptable. This trajectory is shown as a solid 
line in the figure and compared with the DL trajectory, which is plotted as a dashed line. 
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Figure 24: Determination of the Pomeron trajectory from the reaction jp — > (pp. The dots are 
the values of the trajectory at a given —t as determined from Fig. 23 and the full line is the result 



of a linear fit to these values. The Pomeron trajectory as determined by DL [56] is shown for 
comparison as a dashed line. 
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As can be seen from p0| ) and (21), the intercepts and slopes from this determination are in good 
agreement with each other. Whereas the intercepts agree very well with DL, the slopes of the 
present measurement are clearly smaller than that of the DL trajectory, which was determined 
from pp elastic scattering data. It would therefore seem that this simple Pomeron trajectory is not 
universal. 



12 Summary 

Elastic and proton-dissociative photoproduction of p , <ft and J/ifi mesons have been investigated 
at an average photon-proton center-of-mass energy of 94 GeV and for values of —t up to 3 GeV 2 . 
The proton-dissociative event sample was limited to values for the mass of the dissociated proton 
system below 7 GeV. The differential cross section da/dt, for each of the vector mesons in each of 
these two processes, has been measured. The following features are common to these reactions: 

• the ratio of the differential cross sections for the elastic and proton-dissociative reactions 
drops rapidly from a value of about 1 at — t = 0.4 GeV 2 to a value < 0.1 for — t > 1 GeV 2 ; 

• parameterization of the differential cross section as a single exponential, d<r/d£ oc exp(6i), 
for values of —t exceeding 0.5 GeV 2 , yields a difference in the exponential slopes, b, for the 
elastic and proton-dissociative reactions of about 3.5 GeV -2 for each of the vector mesons; 

• the measured ratio of the elastic and proton-dissociative differential cross sections is similar 
to that measured in pp elastic and single-dissociative scattering and is consistent with the 
hypothesis of Regge factorization. 

The analysis of the decay-angle distributions for pion-pair photoproduction in the p° mass region 
indicates a small deviation from SCHC giving non-zero single- and double-flip amplitudes. 

A comparison of the measured differential cross sections, da/dt, for the process 7p — > VN with 
QCD models shows that the perturbative part of the calculations for p° and cf) production [Q 
at the — t values covered in this analysis is well below the data. However, the perturbative QCD 
prediction [11] compares satisfactorily to the J ftp data for values of — t as low as 1 GeV 2 . 



The ratio (j> : p° of the elastic cross sections increases with — t and approaches 2 : 9 at — t ~ 1 GeV 2 . 
The ratio J/ip : p° of the proton-dissociative cross sections increases with — t, but even at —t = 
2.2 GeV 2 is still much lower than the value of 8 : 9 expected for a flavor-independent production 
mechanism. 

The Pomeron trajectory was determined using elastic production of <j> and p° mesons by studying 
the W dependence of da/dt at fixed t values, together with lower W data from other experiments. 
The resulting trajectories are: 

• 7 p -> p°p : a p (t) = (1.096 ± 0.021) + (0.125 ± 0.038)i; 

• ryp^cfrp: a F (t) = (1.081 ± 0.010) + (0.158 ± 0.028)t. 



The values obtained for a p (0) are in good agreement with those of DL [56|. However, the slopes, 
a' , are significantly lower than the value found in pp elastic scattering. 

In conclusion, the results presented in this analysis suggest that even for the highest —t values 
studied here {—t ~ 2-3 GeV 2 ), the variable t cannot be consistently treated as a hard scale in 
perturbative QCD. 
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ZEUS 1995 -f P ^p°N 


-£(GeV 2 ) 


da/dt{pb/GeV 2 ) 


0.353 


8* ±0.733 33 


0.449 


5.i4 ±o.46i8:» ig:» 


0.548 


3.57 ±o.34i818 ±8:S 


0.647 


2.68 ± o.28±8:f? ±8:11 


0.747 


o 97 _|_ n 9R+0.47 +0.52 
Z.Z/ 3= U.ZU_ 068 _ 040 


0.848 


i.62 ± o.2oi8:l! I8:g 


0.949 


i.49 ±o.i9i8:lg ±8^ 


1.093 


i r\n in i r)+0.23 +0.38 

i.ur m u.iz_ 017 _ 027 


1.334 


53 + 06 +al0 +0 - 26 

U.OJ it U.UD_ 07 _ 16 
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26 ± 06+ 010 +0 ' 17 

U.ZD It U.UO_ 05 _ Q]11 


1.937 


o.22 ± o.o4i8:^ 18:58 


2.273 


o.o64 ± o.oi8i8:8g 18:8£ 


2.711 


o.o39 ± o.oioi8:8}l 18:8S 



Table 1: The differential cross sections, da /At , for proton-dissociative p° photoproduction for 
(W) = 94 GeV and M^ < 0.1W 2 . Statistical, systematic and uncertainties due to the modeling 
of the proton-dissociation process are given separately. The normalization error of 15% is not 
included. 



ZEUS 1995 -f P -^p°p 


-t(GeV 2 ) 


da/dt(pb/GeV 2 ) 


0.350 


5.7 ±0.41^ 18:1 
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2 99±0 21 +0 - 62 +0 ' 04 
z.yy it u.zi_ 058 _ 06 
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i.24 ± o.i5i8 j? !8-?I 


0.645 


o.79 ± o.o8i8:l 6 7 18:8^ 
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50 + 06+ 021 + 003 
u.ou ± u.uo_ 10 _ 04 
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23 ± 04+ 013 +003 
u.zo it u.u^_ 05 _ 03 
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0.183 ± 0.033_ : 032 _ : 024 
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o.o53 ± 0.01218:815 18:811 
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Table 2: The differential cross sections, du/di , for elastic p° photoproduction for (W) = 94 GeV. 
Statistical, systematic and uncertainties due to the modeling of the proton-dissociation process are 
given separately. The normalization error of 15% is not included. 
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ZEUS 1995 


-£(GeV 2 ) 


da/dt(>yp->p p) / da/dt(-fp^p°N) 


0.351 
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Table 3: The ratio of the elastic to the proton dissociative cross sections for p° photoproduction for 
(W) = 94 GeV and Mjy < 0.11U 2 . Statistical, systematic and uncertainties due to the modeling of 
the proton-dissociation process are given separately. 



ZEUS 1995 -f P ^p°N 
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Re[r? 4 ] 


r 04 
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Table 4: The spin density matrix elements, rg 4 , Re^ 4 ] and rj 4 in four t intervals for the proton- 
dissociative sample. The data cover the kinematic range 85 < W < 105 GeV and 0.55 < M^ < 
1.2 GeV. Statistical and systematic uncertainties are given separately. 
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Table 5: The spin density matrix elements, r[j 4 , Re[r^ 4 ] and r^ 4 in three t intervals for the elastic 
sample. The data cover the kinematic range 85 < W < 105 GeV and 0.55 < M wn < 1.2 GeV. 
Statistical and systematic uncertainties are given separately. 
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ZEUS 1995 -fp^cpN 
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da/di^b/GeV 2 ) 
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Table 6: The differential cross sections, dcr/dt , for proton-dissociative photoproduction for 
(W) = 94 GeV and Mjy < O.IW 2 . Statistical, systematic and uncertainties due to the modeling 



of the proton-dissociation process are given separately, 
included. 



The normalization error of 15% is not 
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Table 7: The differential cross sections, da/dt , for elastic <fi photoproduction for (W) = 94 GeV. 
Statistical, systematic and uncertainties due to the modeling of the proton-dissociation process are 
given separately. The normalization error of 15% is not included. 
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Table 8: The ratio of the elastic to the proton dissociative cross sections for <j) photoproduction for 
(W) = 94 GeV and M 2 ^ < 0.11U 2 . Statistical, systematic and uncertainties due to the modeling of 
the proton-dissociation process are given separately. 
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ZEUS 1995 jp^J/tpN 
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Table 9: The differential cross sections, da/dt , for proton-dissociative J/ip photoproduction for 
(W) = 94 GeV and M^- < 0.11U 2 . Statistical, systematic and uncertainties due to the modeling 
of the proton-dissociation process are given separately. The normalization error of 15% is not 
included. 



ZEUS 1995 jp^J/ipp 


-t(GeV 2 ) 


dcr/dt^b/GeV 2 ) 


0.219 


o.o96 ± o.o22±g : gs igigg* 


0.567 


0.020 ± o.oo8i8:88^ ±8:8Si 


1.084 


o.ooT ± o.oo5±g:8» ±g : gg 



Table 10: The differential cross sections, da/dt , for elastic J ftp photoproduction for (W) = 94 GeV. 
Statistical, systematic and uncertainties due to the modeling of the proton-dissociation process are 
given separately. The normalization error of 15% is not included. 



ZEUS 1995 


-t(GeV 2 ) 


da/dt(-yp^J/ipp) J da/dti'yp-^J/tpN) 


0.231 


i.3i ± 0.70I8 j° 18:3 


0.596 


i.o8 ± o.97t8:l 18:1 


1.132 


o.i9 ±o.i9i8:3? ±8:g 



Table 11: The ratio of the elastic to the proton dissociative cross sections for J/ip photoproduction 
for (W) = 94 GeV and Mjy < 0.1TU 2 . Statistical, systematic and uncertainties due to the modeling 
of the proton-dissociation process are given separately. 
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ZEUS 1995 -yp^p°N 


-i(GeV 2 ) 


da/dt(fib/GeV 2 ) 


0.352 


4.69 ±0.36 ±0.90 ±0.16 


0.449 


2.61 ±0.23 ±0.42 ±0.15 


0.547 


1.71 ±0.16 ±0.24 ±0.13 


0.647 


1.20 ±0.13 ±0.23 ±0.07 


0.746 


0.96 ±0.11 ±0.27 ±0.05 


0.846 


0.68 ±0.08 ±0.16 ±0.03 


0.949 


0.60 ±0.08 ±0.11 ±0.03 


1.093 


0.39 ±0.04 ±0.08 ±0.03 


1.333 


0.191 ± 0.019 ± 0.036 ± 0.014 


1.638 


0.089 ± 0.020 ± 0.032 ± 0.018 


1.936 


0.071 ± 0.012 ± 0.019 ± 0.010 


2.281 


0.021 ± 0.006 ± 0.008 ± 0.007 


2.714 


0.012 ± 0.003 ± 0.004 ± 0.002 



Table 12: The differential cross sections, da/dt , for proton-dissociative p° photoproduction for 
(W) = 94 GeV and M^ < 0.051U 2 . Statistical, systematic and uncertainties due to the modeling 
of the proton-dissociation process are given separately. The normalization error of 15% is not 
included. 



ZEUS 1995 -yp^cpN 


-i(GeV 2 ) 


da/dt(pb/GeV 2 ) 


0.467 


0.43 ±0.12 ±0.16 ±0.05 


0.639 


0.19 ±0.04 ±0.07 ±0.02 


0.895 


0.096 ± 0.017 ± 0.026 ± 0.004 


1.264 


0.046 ± 0.010 ± 0.014 ± 0.005 


1.977 


0.007 ± 0.002 ± 0.002 ± 0.001 



Table 13: The differential cross sections, da/dt , for proton-dissociative (ft photoproduction for 
(W) = 94 GeV and M^ < 0.051U 2 . Statistical, systematic and uncertainties due to the modeling 
of the proton-dissociation process are given separately. The normalization error of 15% is not 
included. 
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ZEUS 1995 jp^J/ipN 


-i(GeV 2 ) 


da/dt{pb/GeV 2 ) 


0.236 


0.066 ± 0.027 ± 0.033 ± 0.004 


0.616 


0.013 ±0.009 ±0.010 ±0.001 


1.170 


0.023 ± 0.010 ± 0.012 ± 0.001 


2.152 


0.006 ± 0.003 ± 0.003 ± 0.001 



Table 14: The differential cross sections, da/dt , for proton-dissociative J/ip photoproduction for 
(W) = 94 GeV and M^ < 0.05W 2 . Statistical, systematic and uncertainties due to the modeling 
of the proton-dissociation process are given separately. The normalization error of 15% is not 
included. 



ZEUS 1995 


-i(GeV 2 ) 


da/dt(jp^(f>N) J da/dt(~fp^p°N) 


0.492 


0.145 ± 0.028±^ 


1.007 


0.149 ±0.02llg:o2? 


2.047 


0.172 ± 0.056±g:o48 



Table 15: The ratio of the cross sections da/dt for <p to p° for proton-dissociative photoproduction 
for (W) = 94 GeV and M 2 ^ < O.IW 2 . Statistical and systematic uncertainties are given separately. 



ZEUS 1995 


-t(GeV 2 ) 


da/dt('-{p—>-(f)p) J da/dt^p— >p°p) 


0.464 


0.155 ±0.023+^25 


0.865 


0.225 ± 0M7t° o lll 



Table 16: The ratio of the cross sections da/dt for <f> to p° for elastic photoproduction for (W) 
94 GeV. Statistical and systematic uncertainties are given separately. 



ZEUS 1995 


-i(GeV 2 ) 


da/dt(yp^J/ipN) / da/dt(-f P ^p°N) 


0.499 


0.006 ± 0.004±°;°° 2 


1.044 


0.032 ± 0.014±g;g^ 


2.134 


0.097 ± 0.058+g;^ 



Table 17: The ratio of the cross sections da/dt for J/ip to p° for proton-dissociative photoproduction 
for (W) = 94 GeV and M^ < 0.1VU 2 . Statistical and systematic uncertainties are given separately. 
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ZEUS 1995 


-t(GeV 2 ) 


da/dt^p—rJ/t/jp) J da/dt( / -fp—>p°p) 


0.475 


0.015 ± 0.005+g;^| 


0.901 


0.045 ± 0.030±g;gf? 



Table 18: The ratio of the cross sections da/dt for J/ip to p° for elastic photoproduction for 
(W) = 94 GeV. Statistical and systematic uncertainties are given separately. 
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